On the dynamics of a particle on a cone 

K. Kowalski and J. Rembielinski 

Department of Theoretical Physics, University of Lodz, ul. Pomorska 149/153, 90-236 
_ \ Lodz, Poland 

o 

<• 

Abstract 

(N 

^ \ A detailed study of the classical and quantum mechanics of a free particle on 

a double cone and the particle bounded to its tip by the harmonic oscillator 
potential is presented. 



^ 
Ph' 



. , Keywords: Quantum mechanics; Classical mechanics; Constrained 

c^ \ systems; Motion on a cone; Harmonic oscillator; Classical and quantum 

^ . ' unstable motion. 

G^; 

O^ . 1. Introduction 
> 



(N 



O 



S 



The quantum mechanics on conical spaces has attracted much attention 



■^ \ not only in early investigations [ll, |2| but also in recent papers [sl, y, |5|, la, l7| . 



Among different motivations for the study of a quantum particle on a cone 
we only recall the context of the 2 + 1 dimensional quantum gravity [l|, 

cri ■ cosmic strings |8| and defects in various media |9|. Nevertheless, to our best 

knowledge, all investigations of quantum mechanics on a cone were restricted 
to the case of the cone with a single nappe or equivalently the plane with a 

S^ \ deficit angle and the situation of a double cone, i.e. two cones placed apex 

to apex, has not yet been pursued. This is all the more surprising since the 
case of a double cone is much more simple. In particular, in contrast to the 
situation with a single nappe, we have no troublesome questions like that 
concerning behavior of a free particle moving in a generator towards a cone 
tip after reaching it, and there is no need to analyze self-adjoint extensions 
of symmetric operators representing observables of the system. In this work 
we perform a detailed analysis of the classical and quantum mechanics on 
a double cone involving the case of a free particle and harmonic oscillator 
on the cone. An interesting feature of the dynamics on the double cone 
discussed in this paper is the instability of the classical rectilinear motion on 
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the meridian in the sense of Lyapunov and its imphcations on the quantum 
level. The paper is organized as follows. In Sec. II we study the dynamics of 
the free classical particle on the cone. Section III is devoted to the quantum 
mechanics on the cone in the case of the free motion. The classical harmonic 
oscillator on the cone is discussed in Sec. IV. In Sec. V we investigate the 
quantization of the harmonic oscillator on the cone. 

2. Classical mechanics of a free particle on a cone 

Let us consider a particle confined to a surface of a circular double cone 
given by 

xi = / sin a cos (p, 

X2 = Isinasirup, (2.1) 

X3 = I cos a, 

where / G (— cxd, 00) is the coordinate of a particle on a meridian (generator), 
2a is the opening angle of the cone, so a G (0, |), and if G [0,27r) specifies 
the position of a particle on a parallel. In view of (2.1) the equation of the 
conical surface is of the form 

xl + xl-tg^axl = 0. (2.2) 

The Lagrangian of the free particle with mass m constrained to the conical 

surface (2.2) is 

mP msin^aPip'^ ,„ „, 

L = \ —. (2.3) 

2 2 ^ ' 

The corresponding Hamiltonian can be written as 

H = ^ + ^^-^^, (2.4) 

2m 2mP sin a 



dL 



where pi = ^ = ml is the momentum referring to the motion in the merid- 
ian and p<p = ^ = rnl"^ sin^ a^p is the angular momentum. Therefore, the 
Hamilton's equations are 



m' 



^ 



Pi 



Pv 



mP sin^ a' /'2 5^ 

pI 



mP sin^ a ' 
p^ = 0. 



Now let J = Pip = const designates the conserved angular momentum. 
For J = we find from (2.5) if = const, pi = const, and 

l = lo + ^t. (2.6) 

m 

Therefore J = refers to the rectilinear uniform motion along a meridian. 

Consider now the case J ^ 0. Taking into account the following expres- 
sion on the energy 

^ = ^^ + . ,. • 2 ^ (2-7) 



2m/^ sin^ 



a 



we find 



/ = ±^/ — (t + C)2 + ^^-^,-^, (2.^ 



'2E 72 

_(t + C)2 + 
m 2m_E sin^ a ' 

where + (— ) sign refers to the motion on the upper (lower) nappe of the 
cone and C is an integration constant which can be fixed with the help of the 
initial data. Thus, it turns out that for J 7^ the particle resides definitely 
on the upper or lower nappe and there is no communication between these 
two regions of the configuration space. This means, among others, that the 
solution to (2.5) with p^ = J = referring to rectilinear motion across 
the tip of the cone is unstable in the Lyapunov sense. Indeed, an arbitrary 
small perturbation of the initial condition J = leads to large deviations 
of corresponding solutions. Such instability is illustrated in Fig. 1 (bottom- 
right figure). On the other hand, an immediate consequence of (2.8) is the 
following inequality: 

1^1 > ^ , (2.9) 

\'2mE sin a 

which means that for J 7^ the geodesies have the lower (upper) bound on 
the upper (lower) nappe. We remark that the limit J ^ for (2.8) is 

lim/(t) = —t + lo, (2.10) 



J->0 



m 



where we confined for brevity to the case of the motion on the upper nape. 
Clearly this asymptotics does not correspond to any real motion. If so it 
would violate the uniqueness of the initial conditions problem (Cauchy prob- 
lem) related to the system (2.5). An example of the free motion on the cone 
(geodesic) given by (2.1) and (2.5) is presented in Fig 1. Further, from the 



x-i 




x^ 



Figure 1: Top left: the trajectory (geodesies) on the cone (2.1) which is the solution of 

(2.5) with the initial data lo — b, Lp^ = 7r/2, p/o = —1, J ~ ^, and the parameters a = 7r/4 

and m ~ 1. Top right: the rectilinear motion in the meridian of the cone (geodesies) 

referring to J = (see (2.6)). The parameters and remaining initial conditions are the 

same as in the figure on the left. Bottom left: the trajectory given by (2.1) and (2.5) 

with the same parameters and initial data as in the figures above besides of J = 0.01 

corresponding to the limit (2.10) and (2.13). Bottom right: solid line: the solution to 

(2.5) with J = presented in top-right figure. Dashed line: the solution to (2.5) with the 

same parameters and initial conditions as in the top- left figure besides J — 0.1. 
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second equation of (2.5) we get 

2Esma 



and 



arctg- 



J 



/, ^\ 2Esma^ 
(t + C) - arctg C 



ip = const, for J = 0. 

The limit J — )■ for the angle (2.11) is 

V? — (/9o = or \ip — ipol = TT . 

Finally from the first equation of (2.5) and (2.8) we find 

2E{t + C) 



Pi 



2E 



for J ^ 0. 



for J ^ 0, 

(2.11) 

(2.12) 
(2.13) 

(2.14) 



(t + C)2 + 2^^^.^, 



By means of (2.8) pi vanishes when the minimal value is reached by |/| spec- 
ified by (2.9). The limit J — )■ for the momentum in the case of the upper 
nappe can be written as 

lunMt) = \pio\e(t+'^), (2.15) 

where /q > and e{x) is the sign function. 

3. Quantum mechanics of a free particle on a cone 

In quantum mechanics the dynamics of a free particle on a double cone 
defined by the Hamiltonian (2.4), is described by the Schrodinger equation 



.df{i,v;t) pf 



dt 



+ 



J2 



2m 2mP sin^ 



/(^^;^), 



(3.1) 



a 



where pi, I, and J are the momentum, position and angular momentum 
operators, respectively, and we set h = 1. The Hilbert space for a quantum 
particle on a cone is specified by the scalar product 



{f\9)= / d^ dl\l\ni,^)gil,^), 

Jo J~oo 



(3.2) 



where \l\dldip coincides up to the multiphcative constant with the surface ele- 
ment of the cone defined by (2.1) such that dS = sma\l\dld(p. The operators 
pi, I and J act in the representation (3.2) in the following way: 

if{l,v) = lf{l,^), (3.3) 

We remark that the action of the operator pi different from the Schrodinger 
representation ensures its hermicity with respect to the scalar product (3.2). 
Of course, such form of the operator pi preserves the canonical commutation 
relations of pi and I as well. The counterpart of the operator pi in the case 
of the cone with a single nappe was introduced in reference 4. Nevertheless, 
it is not self-adjoint so it cannot be regarded as a physical observable. We 
stress that in opposition to the case of the single nappe (or the plane with 
deficit angle) where the observables such as for example energy are labelled 
by parameters related to their self-adjoint extensions j^, l5|, the operators 
representing physical observables for the quantum mechanics on the double 
cone discussed herein, are defined uniquely. We finally point out that the 
analysis of self-adjoint extensions in the case of the cone with a single nappe 
is closely related to the study of self-adjoint extensions describing a quantum 
particle on a plane with extracted point performed by us in reference 10. 

Consider now the eigenvalue equation for the Hamiltonian Hf{l, ip) = Ef{l, ip) 
following directly from (3.1) and (3.3) 

1 / a^ 1 (9 1 1 9^ , 

' fE{l,^) = EMl,^). (3.4) 



2m V dP I dl 4/2 p sin^ a dip^ 

On separating variables we find 

/,,,,(/,¥.) =e^^%-s(/), (3.5) 

where fj^Ei^f) are the common eigenvectors of the Hamiltonian and the 
operator of the angular momentum J, and Uj^e{1) satisfies the equation 

-^P- + 7^^ + 2mE - -^^ u,,e{1) - 0. (3.6) 



The solution to (3.6) for j 7^ with convergent norm (in a distributive sense) 
is of the form 

uj,e{1) = CJ r^-^{V2^l), (3.7) 

V 4 + sin^a 

where Ji,{z) is the Bessel function of the first kind and C is a normahzation 
constant. Now, bearing in mind the behavior of the classical free particle for 
J 7^ discussed in previous section, taking into account that Ju{0) = for 
z/ > 0, and the identity 

-M-z) = e^^'Mz), (3.8) 

it seems plausible to identify the case / > (/ < 0) in (3.7) with the quan- 
tum particle on the upper (lower) nappe. Let us designate by m+j,e(/, ip) 
{u^j^Ei},^)) the corresponding solution to (3.6), so 



-(V2rnEl), for / > 0, 
u+je{1) = { Vt+sii^ (3.9) 

for / < 0, 




for I > 0, 



{y/2mEl), for / < 0. 



In general case the state of a particle on a cone is the superposition of the 
states M+j,_B and M-j,_b, so we finally arrive at the solution to (3.6) such that 

/,, s(/, ^) = e'^^[Au+jM + Bu^,,e{1)], J + 0, (3.11) 

where A and B are constant. Taking into account the so called "closure 
equation'' 



11 of the form 



1 1 

xJ^{ax)Jy{(5x)dx = —6{a — (3), v > — , (3-12) 

a 2 

and (3.8) we find for the eigenvectors of the Hamiltonian (3.11) normalized 
as 

/'27r roo 

ifjM'E') = d^ dl |/|/^^(/, ^)fj>,E'{l. ^) = Sj,,SiVE - v^), 

Jo J-00 

(3.13) 
the relation 



, .,9 , ,9 my E 

|A|2 + |fi|2 = ^_. (3.14) 

TT 



Now, it can be easily checked that for j = the eigenvectors of the 
Hamiltonian satisfying (3.6) are given by 

fo^E{l,^) = ^=[Asm{V2mEl) + B cos{V2mEl)]. (3.15) 



For E"^ + E'"^ 7^ the coefficients A and B of the normahzed state (3.15) 
satisfy the condition jAp + \B\'^ = \/2m/27r^. Otherwise, ii E = E' = we 
have \B\ = ^. Taking into accout the identity 

[Y 

Ji{z) = \ — sinz, (3.16) 

we find that in the hmit j — ?■ of (3.11) only a part of the solution (3.15) 
referring to the rectilinear motion in the meridian is obtained. The lacking 
part could be derived as the limit j — )■ of the function J n — -2—{V2mEl) 

by means of the identity 

J-i{z) = \ — C0S2;. (3-17) 

2 \ TIZ 



However, the Bessel function J n — -2— (v2mi^/) has the index lesser than 

~V4+— 2-^ 

— I and the "closure equation" (3.12) related to normalization cannot be 
applied. Such behavior of solutions (3.11) and (3.15) is a quantum scar of 
the instability of the classical rectilinear motion in the cone discussed in Sec. 
II. 



4. Classical mechanics of the harmonic oscillator on a cone 

We now discuss the motion of the classical particle on the double cone 
bound to its tip by harmonic oscillator potential. The coresponding Hamil- 
tonian is (see (2.4)) 



H^± + —^ + '^e, (41) 



and the Hamilton's equations are of the form 

Pi_ 
m 



I 



V 



Pi 



Pv 



mP sin^ a ' 



(4.2) 



mP sin^ a 
p^ = 0. 



— mu I, 



Of course J = p^ = which leads to ip = const, refers to the case of the 
standard harmonic oscillator on a line (generator) with the solution 



/ = Iq cos ut 



Pio 
mu 



sin cut, 



Pi = Pio cos ut — uitlIq sin ut. 
Now, using the expression on the energy such that 



(4.3) 



mP 

E=-- + 



J2 



2mP sin^ a 



+ 



mu 



I' 



(4.4) 



we find for J 7^ 



l = ± 



\ 



E 



E^ - ^sin2u{t + C) 



mu" 



J^O, 



(4.5) 



where plus and minus sign correspond to the upper and lower nappe, respec- 
tively, and C is an integration constant. It thus appears that as with the 
case of the free motion, whenever J 7^ the particle remains in one concrete 
nappe during the time evolution. Therefore, the solution (4.3) describing the 
harmonic oscillations in a meridian around the point / = turns out to be 
unstable (see Fig. 3). On the other hand, an immediate consequence of (4.5) 
is that for J 7^ the motion is finite. Namely, we find 



\ 



E 



E^ 



J2w2 



< 



< 



mu" 



\ 



E 



E^ 



J2cj2 



mu" 



(4.6) 



We remark that in the limit J — )■ we obtain the correct amplitude \ h^^ of 
the harmonic oscillations (4.3), however the period of harmonic oscillations 
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(4.3) with J = is twice the period of oscillations (4.5) with J 7^ 0. For the 
upper nappe the limit J — )> of (4.5) is 



The trajectories on the cone (2.1) satisfying (4.2) involving the case of small 
J are shown in Fig 2. 

Furthermore, using the second equation of (4.2) we get for J 7^ 



e{J) 


/ Eiguj{t + C) + J E^ - 

arrtp- 




EtgcoC + Je^ - 

- arrtcr 




sma 


1 "" LPu.^ 


/ J2t.2 
V sin*^ a 





\ y sin^ a y sin 

(4.8) 
As with the case of the free motion we obtain in the limit J — )■ 

ip — ipo = or \ip — ipol = n. (4.9) 

Finally, using the first equation of (4.2) we derive the following formula on 
the momentum 

Pi = ± : ' ^'"'" ^ • ' . (4.10) 



+ 


^M^JE^- 


sin*^ a 


cos2a;(t + 


C) 




\h + V 












Ie- 


.Pu. 
sin*^ 


-sin2w(t 


+ C) 



We point out that pi is zero when the extremal values of |/| are reached given 
by (4.6). In the limit J — )■ this expression takes the form 



fo^^cos2a;t-f^-l)sin2a;t 



Pi = +^A^^ — ' ^ ^ (4.11) 

^fii sin 2oot+('^-l) cos 2ujt 



5. Quantum harmonic oscillator on a cone 

Consider now a quantum particle moving on a double cone and bound to 
its tip by the harmonic oscillator potential. The corresponding Hamiltonian 
is 

H = ^-^ + + 1^. (5.1) 

2m 2mPsm^a 2 
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x, 



X, 




Figure 2: Top left: the trajectory of the harmonic oscillator on the cone given by (2.1) and 
(4.2) subject to the initial data Iq = 9, (/?o = 0.1, piQ — —1, J = 20, and the parameters 
a = IT I A, m = 1 and ui = \/2- View from below. Top right: view of the right side of the 
trajectory on the cone presented in the figure on the left. Bottom left: the motion of the 
harmonic oscillator on the cone with the parameters and initial conditions the same as in 
the top figures besides of J = 0.01 referring to the limit (4.7) and (4.9). View from below. 
Bottom right: right-hand view of the bottom- left figure. 
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xj - 




Xi - 




Figure 3: The illustration of the instability of the harmonic oscillations on the cone cor- 
responding to the vanishing angular momentum. Left: white line refers to the solution to 
(4.2) with the same values of parameters and initial conditions as in the top figures in Fig. 
2 besides of J = 4. Black line is the part of the segment from the lower nappe which is 
the trajectory of the harmonic oscillator with J = (the remaining parameters and initial 
data coincide with those for the white trajectory). Right: the view of the figure on the 
left from the behind. 
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Consider the eigenvalue equation for the Hamiltonian 

HfE{l,ip) = EfEil,^). (5.2) 

On separating variables and setting 

/,,^(/,(^) = 4^e'^%.^(0, (5.3) 

we can write (5.2) in the form 

^%^ + (2mE - -l^ - [mufi] uMl) = 0. (5.4) 

at'' \ [^ sm a J 

We remark that the potential occuring in the Schrodinger equation (5.4) 
which is a sum of the harmonic oscillator and inverse square potential is 



referred to as the Smorodinsky-Winternitz potential [12|, isotonic oscillator 



potential [13[ or harmonic oscillator with "centripetal barrier" vM- This 
potential is also closely related to the well-known Morse potential jl5|. How- 
ever, as far as we are aware, the Schrodinger equation (5.4) was not discussed 
in the literature in the context of the double cone. The quantum harmonic 
oscillator in the case of a cone with a single nappe or equivalently the plane 
with a deficit angle was studied in jj] and |6|. As mentioned earlier this case 
is much more complicated than that investigated herein dealing with quan- 
tum mechanics on a double cone, because one is forced to analyse self-adioint 
extensions of symmetric operators representing physical observables \4, l5| . 
Now, making the ansatz 

Uj,e{1) = C|/|'*e~^^w(mcu/^), (5.5) 

where C is the normalization constant, we find that whenever the parameter 
s satisfies 

s(s-l) = -^, (5.6) 

sm a 



which leads to 




(5.7) 



then w fulfils the Kummer equation 



d'^w ( 1 \ dw f s 1 £■,,, ^ ,^, 

13 



A numerically satisfactory solution to (5.8) with a good behavior at infinity 
is given by the Kummer function 16 1 



w{x) = U[- + ---,s + -,x 



(5.9) 



Another notation for the Kummer function U{a,b,z) is \E'(a, 6, z) 17|, Il8 
Now, U{a,b,x) is a polynomial when a = —n, where n = 0, 1, 2 . . ., and 
U{a, b, x) ~ x~°- for a; — t- oo. Demanding decreasing of Uj^Ei}) as / — )■ oo this 
leads to the quantization condition on energy 



E, 



3,ri 



s 1 

2uj{n + - + - 



(5.10) 



Furthermore, the requirements J 7^ and convergence of the norm of the 
solution to (5.8) rule out the minus sign in (5.7), so 



E. 



],n 



1 1 I 4?2 



sin^a 



(5.11) 



and 

w{x) = U {~n,l + ^\ll 

Using the identity [1 



4j^ 



2, V sm a 



X 



U{-n,a + 1, z) = {-l)^n\L^{z), 
where L^{x) are the generalized Laguerre polynomials, we get 



wix) 



-irniu 



'1+ 



4jZ 



(5.12) 



(5.13) 



(5.14) 



Hence, we finally obtain the desired normalized (in the sense of L^(]R, da;)) 
solution Uj^nil) to (5.4) such that 



Mj,n(0 = ^i,Sn(0 



-IT 



\ 



(muoY^y^^^nl 



i(iWi+-^ 



r n+1 



1 + 



4j2 



X L 



Wi+;S 



'" "(mwr). 



(5.15) 
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and the normalized solution to the Schrodinger equation on the cone (5.2) of 
the form 

/,>(/, ^) = -i— e'^%- „(/), J ^ 0, (5.16) 

lTl\/\l\ 



where the normalization is given by 

/»27r /'oo 

{fj,n\fj'n')= dip dl\l\fl^{l,ip)fj'y{l,ip)=5jf5nn'- (5.17) 

Jo J -oo 

Up to the phase factor and normalization constant the solution (5.15) coin- 
cides with the solution to (5.4) introduced by Hakobyan et al in reference 12. 
However, the parameter k labeling the solution obtained in 12| had no phys- 
ical interpretation. On the other hand, it follows immediately from (5.15) 
that if /c > |, then this parameter can be related to the angular momentum 



7 on the cone via k = 7t\/ 1 + . ^ • 

Now, an immediate consequence of (5.3) and (5.4) is the following nor- 
malized solution to the eigenvalue equation (5.2) for j = 0: 





/o, 


{]\ — 


1 


n 
V 


TT 


'-r 


1 


-.e 


mujl 




'" 2.V1I1 


V2"n 




referring 


to the 


energy 

















HniVmuJl), (5.18) 



Eo,n = u;{n + i), (5.19) 

where Hn{z) are Hermite polynomials. Indeed, (5.4) for j = is the Schrodinger 
equation for the standard harmonic oscillator on the line. Analogously as 
with the free motion, the solution (5.18) corresponding to j = is different 
from the limit j — ?■ of the solution (5.16). In fact, taking into account the 
identity [19 

H2n+i{x) = (-l)"22"+ir2!a;L|(x2), (5.20) 

we obtain the limit j — t- of fj^n{l) given by (5.16) such that 

V f n\ 1 (muj\\ 1 _i2i^„ , , ,,,. 

,•^0^^'"^^ 27rv^V7ry ^22"+i(2n + 1)! ^^ ' '^ 

(5.21) 
The energy spectrum in the limit j — )■ is 

limE,„ = o;(2n+l + |). (5.22) 

j->-0 
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We have thus obtained in the hmit J — )■ only a part of the solution cor- 
responding to the case of the harmonic oscillator on a meridian (generator) 
related to the odd Hermite polynomials. As in the case of the free motion 
this observation is consistent with instablility of the solution (4.3) of classical 
equations of motion describing harmonic oscillations around a point / = in 
the meridian. 

6. Conclusion 

In this work we study the classical and quantum particle on the double 
cone in the case of the free motion and harmonic oscillator potential. An 
interesting pecularity of the dynamics on the cone is the instability of the 
free motion in the generator on both classical and quantum level. This ob- 
servation seems to be of importance for better understanding of quantum 
dynamics. Indeed, as with deterministic chaos the role of instabilities in 
quantum dynamics is far from clear. Referring to chaotic motion it is also 
worthwhile to point out that the unstable solution describing harmonic os- 
cillations around a tip of the cone is rather exotic example of the unstable 
and bounded solution of a nonlinear dynamical system which does not show 
chaotic behavior. Furthermore, a remarkable property of the dynamics on 
the double cone investigated in this paper is the role played by the vertex 
which acts as a filter selecting from all possible motions the rectilinear one 
in the generator. Finally, we remark that the results of this paper provide 
a physical interpretation as a harmonic oscillator on the double cone for the 
one-dimensional Smorodinsky-Winternitz potential which can be connected 
with the Morse potential. 
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